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Abstract 

We provide the analytic predictions for the flavor non-singlet scalar correlator, which 
will enable determination of the scalar meson mass from the lattice scalar correlator. We 
consider simulations with 2+1 staggered sea quarks and staggered or chiral valence quarks. 
At small u/d masses the correlator is dominated by the bubble contribution, which is the 
intermediate state with two pseudoscalar mesons. We determine the bubble contribution 
within Staggered and Mixed Chiral Perturbation Theory. Its effective mass is smaller than 
the mass of irr/, which is the lightest intermediate state in proper 2+1 QCD. The unphysical 
effective mass is a consequence of the taste breaking that makes possible the intermediate 
state with mass 2M n . We find that the scalar correlator can be negative in the simulations 
with mixed quark actions if the sea and valence quark masses are tuned by matching the 
pion masses M va i >va i = M n5 . 

1 Introduction 

The Nature of the lightest observed scalar resonances is still not revealed yet. In this paper 
we consider the flavor non-singlet scalar meson and the problems related to its simulation 
on the lattice. Below 2 GeV there are two experimentally well established scalar resonances 
with isospin 1=1: ao(980) and ao(1450) [Q. It is still not clear which of the two corresponds 
to the lightest qq scalar meson. Most of the models [2] and lattice simulations have difficulty 
relating the mass and decay properties of ao(980) with the qq state. If ao(980) belongs to 
the lightest qq nonet, it does not have an obvious strange partner since the mass of iTo(1430) 
is relatively high, while k(800) is experimentally very controversial PjJ E] and probably too 
light to be the partner. This suggests the possibility that a (1450) might be the lightest qq 
state with 1 = 1, while ao(980) could be exotic state such as tetraquark qqqq [I] or mesonic 
molecule. 
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The issue could be settled if the mass of the lightest qq state with 1=1 would be reliably 
determined on the lattice. We refer to this resonance as a and we consider du scalar meson 
for concreteness. In order to determine a$ mass, lattice simulations evaluate the connected 
scalar correlator 

C{t) = J2(0\d(x,t)u(x,t) u(0,0)d(0,0)|0) . (1) 

X 

If cto is the lightest state of the Hamiltonian with quantum numbers J p = + and 1 = 1, then 
C(t) oc e~ ma0t at large t and extraction of m a0 is straightforward. Multi-hadron states with 
J p = + and 1=1 also propagate between the source and the sink of the scalar correlator in 
Fig. ^and they may shadow the interesting part e~ ma0t . The most important multi-hadron 
state is the intermediate state with two pseudoscalars P1P2 , we call it the bubble contribution 
B and display it in Fig. 121 The state P1P2 has discrete energies E n on the lattice with spatial 
extent L 



M Pl + M P2J y M Pi + |f| 2 + y Mp 2 + |f | 2 , y M Pi + \ 2 fn\* + + |f n\* (2) 

and gives rise to a sum of decaying exponentials e~ Ent in the correlation function. If the 
masses of P\ and P2 are small, the bubble B presents a sizable contribution to the scalar 
correlator (Fig. ^) and it has to be incorporated in the fit of the lattice correlator (0) 

C(t) = Ae- ma0t + B{t) + ■■■ . (3) 

The dots represent contributions of excited physical scalar meson and other multi-hadron 
intermediate states, which are both less important. 
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Figure 1: The scalar correlator with 1 = 1 receives the contribution from the propagation of the 
physical scalar meson ao as well the propagation of multi-hadron states. The most important multi- 
hadron state is the two-pseudoscalar intermediate state, which gives rise to the so-called bubble 
contribution. 

Let us review the importance of two-pseudoscalar intermediate state in theories with 
various numbers of dynamical quarks: 

• The interest of this paper is 2 + 1 flavor QCD (m u = m ( ij^m s ), where the intermediate 
states are 7177, KK and irn' (inr is not allowed by Bose symmetry and conservation of J p 
and I G in the strong decay). In Nature, all these intermediate states are lighter than 
cio(1450), whereas 7177 is also lighter than ao(980). Hence, two-pseudoscalar intermediate 
states dominate correlation function at large t if simulations are performed at small 
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u/d masses. This makes it difficult to extract the mass of ao from (jSJ). For the same 
reason, the scalar correlator is very sensitive to the possible unphysical approximations 
undertaken in the simulations, like the application of the staggered fermions or mixed 
quark actions. In the present paper we calculate the bubble contribution for the case of 
simulations with staggered sea and staggered valence quarks, as well as for simulations 
with staggered sea and chiral valence quarks. The bubble contribution incorporates 
the physical intermediate states irr], KK, ttt]'. It also gives the size of the unphysical 
effects of staggered fermions and mixed quark actions in the scalar correlator. We 
find that these unphysical effects are sizable. Our results are relevant for the existing 
simulations of scalar mesons with staggered sea and valence quarks by MILC [S] and 
UKQCD jHj. They are also relevant for the simulations with staggered sea quarks and 
chiral valence quarks, currently performed by LHPC, NPLQCD and UKQCD 

• Two flavor QCD allows only irrj' intermediate state, where r( is two-flavor singlet with 
mass M^, = M% + \rn^. The nr/' state is relatively heavy and does not shadow so 
significantly the interesting part of the correlation function given by e~ ma0t . The mass 
of ao was extracted from the conventional exponential fit in the simulation with two 
flavors of dynamical Domain Wall Fermions The resulting mass was unaffected if 
the irrj' intermediate state was also taken into account in the analysis of the correlation 
function jHj. The SCALAR [Oj and UKQCD ^H] Collaborations also simulated ao 
in two-flavor QCD. Partially quenched scalar correlator with two dynamical Domain 
Wall Fermions was simulated in [Hj. The prediction of the bubble contribution within 
Partially Quenched ChPT describes well the striking effect of partial quenching, and 
the scalar meson mass was extracted by fitting the partially quenched correlators to 

• The scalar correlator in quenched QCD is dominated by irrj' intermediate state at light 
quark masses This has large unphysical effect which can be attributed to the 
breakdown of unitarity in quenched QCD. The scalar correlator is negative and has 
the effective mass 2M n at large t. The scalar meson mass m a o has been determined by 
fitting the quenched scalar correlator to the sum of ao exchange and 7tt/ exchange 1 ((HJ) 

in mug 03. 

So the bubble contribution (Fig. EJ) to the scalar correlator is sizable and has to be 
incorporated in the fit of the lattice correlator at light quark masses. Our analytic prediction 
for the bubble contribution B(t) offers a way to determine the scalar meson mass m a0 by 
fitting the lattice scalar correlator to Eq. Q. The extraction of m a0 in this way is possible 
for the range of quark masses and times, where e~ ma0t is not negligible with respect to the 

1 The quenched result of ^3] agrees with two-flavor dynamical result [5j if the effect of quenching is 
incorporated at the leading order in the chiral expansion. The quenched m a o from |13| is somewhat lower 
since the quenching effect is incorporated at the next-to-leading order, given by the diagrams in Fig. 8 of 
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Figure 2: The bubble contribution to the scalar correlator. Here P\ and P2 denote pseudoscalar 
fields in the relevant version of ChPT. 



bubble contribution 2 . We address simulations with 2+1 staggered sea quarks and staggered 
or chiral valence quarks: 

• In the case of staggered valence quarks, we determine the bubble contribution using 
Staggered Chiral Perturbation Theory (SChPT) [14] . All the parameters entering 
our prediction, have been determined by MILC simulations. We find positive bubble 
contribution with effective mass 2M n5 < m e ff < M n5 + M V5 at large t. This offers a 
natural explanation why the MILC j^j and UKQCD [Bj Collaborations observed the 
effective mass well below M n5 + M V5 , which would be the lightest mass in proper 2+1 
QCD. We propose that this unphysical result is due to the taste breaking. 

• In case of chiral (Ginsparg- Wilson) valence quarks we apply Mixed Chiral Perturba- 
tion Theory (MChPT) ^B]- The result is relevant for simulations of LHPC, NPLQCD 
and UKQCD [7] which employ staggered sea and chiral valence quarks. The bubble 
contribution depends on one unknown parameter a 2 AMi X , which gives magnitude of 
the taste breaking in the mass of the pion composed of one sea and one valence quark: 
Mvaisea = ^o( m vai + m sea) + 0? A Mix . We find that the use of the mixed quark ac- 
tions together with the taste breaking can have sizable unphysical effects on the scalar 
correlator and they can make the scalar correlator negative. Namely, the bubble con- 
tribution is found to be negative for c^Amix ^ —0.01 if sea and valence quark masses 
are tuned by imposing M va i tVa i = M JT5 (applied by LHPC Collaboration [7j). The bub- 
ble contribution is positive for the reasonable values of a A A.Mi X if M va i }Va i = M 7TI is 
imposed. 

The bubble contribution to the scalar correlator is calculated in Section 2. It is evaluated 
using Staggered and Mixed ChPT in subsection 2.1 and 2.2, respectively. The implications of 
our analytical results for the lattice simulations are also discussed in this section. Conclusions 
are given in Section 3. 

2 The bubble contribution completely dominates the correlator for very small quark masses and large 
times. 
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2 Contribution of two-pseudoscalar intermediate state 



In this section we calculate the contribution of two-pseudoscalar intermediate state to the 
scalar correlator with a point source and a point sink (0). Our calculation is carried out at 
the lowest order of the effective field theory and the relevant bubble diagram is shown in Fig. 
121 For this purpose we only need the coupling of the scalar current to a pair of pseudoscalar 
mesons P1P2 and the propagators for pseudoscalar s Pi and P 2 . 

The point scalar current is expressed in terms of two pseudoscalar fields $ using the 
Chiral Perturbation Theory ^T] 3 in the Appendix A: 

d(x)u(x) ~ B [$(x) 2 U (4) 

So the coupling of the point scalar current to two pseudoscalars is equal to the slope param- 
eter Bq = M%/(2m q ) which can be determined from the lattice data. The same coupling 
applies also in the Staggered as well as in the Mixed version of ChPT at the lowest order. 
The bubble contribution in Fig. |2]is equal to 

B = (0\du ud|0W = £ 2 (0|[$ 2 U$ 2 W|0> . (5) 

The Wick contractions lead to the products of two propagators for pseudoscalar fields. The 
propagators and the resulting bubble contribution will be given separately for Mixed and 
Staggered ChPT in two subsections that follow. 

We note that the bubble contribution (J3J) applies for the bare lattice point-point correlator 
(JTJ via (jHJ), if the bare quark mass is inserted to B = M%/(2m q ). This is possible since the 
scalar current and Bq oc 1/ m q are multiplied by the same factor Z$ in the renormalization. 

2.1 Staggered sea quarks and staggered valence quarks 

First we consider the simulations with staggered valence quarks and three flavors of staggered 
sea quarks with m u = ^ m s . We study only the dynamical correlators, where the valence- 
quark mass is equal to the u/d sea-quark mass. Staggered quarks have four tastes and the 
scalar correlator is conventionally calculated using the taste-singlet source and sink. This 
correlator receives contribution from the physical scalar meson and the bubble contribution 
(J3J). It also receives an oscillating contribution from a pseudoscalar meson 4 , but the oscillating 
contribution will not be considered here since it has been identified and removed in the 
analysis of the lattice data [5J E] • 

The bubble contribution (0) in a theory with four tastes and three flavors of sea quarks 
(4 + 4 + 4) is equal to 

Bf^I = \Bl <°I[ $2 W «*|0> , (6) 

t,t'=i 

3 B was denoted by r in and by 2/x in j^j. 

4 The operator with spin / and taste I couples also to state with spin 7574 and taste 7574. 
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where 1/4 is appropriate normalisation. Here $ is 12 x 12 pseudoscalar matrix of Stag- 
gered ChPT [T3j and the subscripts u,d denote its flavor component, whereas t,t' = 1, ..,4 
denote its taste component. The field $ = 5Z&=i \T b <p h is expressed in terms of mass 
eigenstates <fr b , where <p h is 3 x 3 pseudoscalar matrix with flavor components <p b j; and 
T b = {£5, i^s^, i^&z, £/}. The bubble contribution can be expressed in terms of the 
pseudoscalar propagators (4> b \4> b ) by applying the Wick contractions to (JHJ) 

Bigg = l^o 5^ (<0LkL> + <0L10L> [<0LI0L> + <&l&> +2<0Ll&>]) • (7) 

6=1 ^ ' 




(a) d (b) d (c) 

Figure 3: The quark level diagrams that correspond to the bubble contribution (Fig. |2j) in 
simulations with dynamical u, d and s quarks. The diagram (c) represents also the diagrams with 
any number of quark loops inserted in the place of the dots. 



The propagators (<f) b \(f) b ) for the pseudoscalar fields of various tastes b are provided by 
SChPT jTJj and are collected in the Appendix B. Propagators for all tastes (J, V, A, T, P) 
have connected part, while tastes /, V and A have also disconnected part. The bubble 
contribution J7J) can be represented by the quark-level Feynman diagrams in Fig. |3] and it 
is calculated by inserting the mesonic propagators in Appendix B to expression (0). We get 
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in 2 • 4 + 4 theory 
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m „w(t a 2 5' VA ) in 2 + 1 theory 



(9) 
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and the functions m v , are given by (J2*3*J) in Appendix B. The four tastes per sea quark 
(2-4 + 4) are reduced to one taste per sea quark (2 + 1) by applying \/Det trick in staggered 
lattice simulations (5J. In SChPT ^3] this is achieved by multiplying every sea quark loop 
by factor 1/4. The sea quark loops are present in diagrams (a,b) and in hairpin vertex of 
diagram (c) in Fig. El The corresponding factors of 1/4 are incorporated in Q. 

The pseudoscalar masses M Ub = M nb = 2B m u + a 2 A b , M Sb = 2B m s + a 2 A b and M USb = 
Bo(m u + m s ) + a 2 A b of unmixed mesons have been determined by MILC simulations JSJ. The 
taste breaking a 2 A b is flavor independent and vanishes in the continuum limit. The hairpin 
couplings of the taste-vector and taste-axial mesons were also determined by MILC for their 
coarse and fine lattices ^3] 

(a 4 ^) coarse = -0.04(1) , (a%) coarse = -0.01(3) , (a 2 c^ y ) finc ~ 0.35 (a 2 5' AV ) cowse . 

(10) 

These hairpin couplings mix the flavor neutral mesons of tastes V and A to the mass eigen- 
states rjy,A and r]' VA with masses given in (|9|23j) . 

a) Results in case of 4 tastes per quark 

The theory with four tastes per sea quark (i.e. without \/Det trick) is not relevant to 
QCD. In spite of that, it is illuminating to consider the properties of the bubble contribution 
in this case before proceeding to the single-taste theory. 

The four-taste theory is a proper unitary quantum field theory even at finite lattice 
spacing and B^+^it) (jBJ) is positive definite for any values of the input parameters. This 
is explicitly demonstrated for the simplified case of m u = = m s and A& = with general 
values of hairpin couplings a 2 5' VA in Appendix C. 

The I?!2+^ T (t) oc e~ 2Mnt at large time for zero and non-zero lattice spacing. This is in 
sharp contrast to three flavor QCD, where B(t) oc e -W*+M v )t with M 2 = j_ M 2 + | M |_ It ig 
easy to understand why an intermediate state with mass 2M n is possible in -Bf^+f 5 " (jBJ) at 
a = 0. The mass eigenstates that enter as intermediate states are ir b rj b and K b K b . Only the 



mass of rji is lifted by the anomaly to y|M 2 + §M|, while the etas of other tastes retain 
mass M n at a = 0. So the lightest intermediate state with 1 = 1 has the mass 2M n . 

b) Results in case of one taste per sea quark 

Finally we study the properties of B^z^ 1 ^ (JBJ), which is relevant to QCD since the four 
tastes of sea quarks have been reduced to one via (JUJ). 



2 A b —> 0, ' 

reduces to the contributions of iri] and KK from 2 + 1 ChPT 



In the continuum limit there is no taste breaking (a A b — > 0, a 5' VA — > 0) and 



T>ChPTf \ _ n2V^ J 2 1 1 1 1 

KP) 0Z^ 3 (fc + p )2 + M 2 k 2 + M 2 T (fc + p y + M 2 k k 2 + M 2 



tj2 /o -(M n +Mn)t -2M K t\ 

lim B c hPT{t) = B± ( + £_ , (11) 



t-oo w 4L 3 V3 M^M^ M 2 K 
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where M 2 = |M 2 + §M§ and M 2 K = |M 2 + |M| The bubble contribution is positive and 
drops as e -( M ^+ M v) t a t large t, as expected in QCD. 

The fourth-root trick destroys the unitarity of the theory at finite lattice spacing and 
^2+i PT (^) can be negative for specific values of input parameters 5 . On top of that, the 
taste breaking at nonzero a allows for the intermediate state with 2M T , which is prohibited 
in proper QCD. This can be understood by considering the quark flow diagrams in Fig. El 
The diagram in Fig. exhibits the expected time-dependence e _2Mx * in ChPT as well as 
in SChPT. The diagrams in Figs. an d Efc both contain terms with e~ 2Mnt , but these 
exactly cancel in ChPT with no taste breaking. This cancellation does not occur in the 
presence of taste breaking when pions of all tastes ^t,i,v,a flow in diagram Eb, but only the 
pions t^i, v,a flow in the diagram Efc. The taste breaking therefore significantly modifies the 
time-dependence 

r^srhPT, s Bl ( 2 e~ 2M ^ 1 terms with pairs \ 

1 m B^ hPT (t) = — 2- — + 2M + F - (12) 

i-oo 2+1 w 4L 3 V16 M 2 5 e ~ 2M ^ tit], irr)', KK J K 1 

with respect to the expected behaviour in proper QCD (jlljl . The coupling to the un-physical 
intermediate state with mass 2M n vanishes only in the continuum limit. 

Now we evaluate B2+i PT (0) for the case of MILC configurations in order to see if our 
observations agree with the results from the lattice simulations. Fig. 0] shows the magnitude 
of the bubble contribution evaluated for the 2 + 1 MILC coarse and fine lattices, using the 
pseudoscalar masses, the bare quark masses, the hairpins (fTUJl and the lattice volumes from 
0. The Euclidean loop momentum k is summed over the discrete momenta in the finite box 
of the lattice 6 . The bubble contribution is larger for lighter u/d quark masses. It becomes 
the dominant contribution to the scalar correlator at large t if two-pseudoscalar state is 
lighter than ao, which is realized 7 for the u/d quark masses considered in Figs. EJandEJ 

The effective mass of B2+i PT {t) is shown by triangles in Fig. It is compared with 
M f5 + M V5 , 2M 7r5 and with the effective mass of the scalar correlator from MILC simulation 
j^]. The effective mass of the analytic prediction lies between 2M 7r5 < m e ff < M ns + M m 
at large t. This agrees 8 with the effective mass observed in lattice simulations by MILC 
(shown in Fig. EJ) and by UKQCD j^J. The exception are the highest quark masses, 
where the lattice simulations |H] observe effective mass even below 2M 7rg , which can not 
be understood within given theoretical framework. The unphysical effective mass 2M n is 
lifted in the continuum limit to the physical value M n + M v , which is the effective mass of 
B ChPT (jllj) and is represented by circles in Fig. |3J 

5 For example, -Bf+ 1 ' lPT (t) is negative in the case of input parameters relevant for MILC simulations with 
sign of a 2 S' v A in fTTijl reversed. 

6 The detailed expression for the finite sum in B(t) is given in Appendix B of [H| for analogous case of 
Partially Quenched theory. 

7 This applies for observed m a o = 1450 MeV as well as m a o obtained by MILC from heavier quark masses. 

8 Fig. presents the effective mass of Sf+' 1 ' lPT at t = 20 for coarse and at t = 30 for fine lattice, while 
the MILC effective mass was extracted from a wider time-range t ~ 5 — 20 [Sj, so they are not expected to 
agree exactly. 
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coarse MILC simulation fine MILC simulation 




Figure 4: The bubble contribution B2?x PT (t) © for MILC simulation with 2 + 1 staggered sea 
quarks and staggered valence quarks. The am u / ( i/am s denote the 2 + 1 sea-quark masses in lattice 
units. The coarse and fine lattices have the volumes 20 3 x 64 and 28 3 x 96, respectively. The thiner 
lines denote the variation of B2+i PT (t) due to the uncertainty in the hairpin parameters ()1UI) . 



2.2 Staggered sea quarks and chiral valence quarks 

The simulations with chiral valence quarks on the available staggered MILC configurations 
present an appealing possibility. The LHPC and NPLQCD Collaborations use domain wall 
valence quarks, while UKQCD uses overlap valence quarks [7]. The corresponding effective 
theory has been derived recently ^H] and we refer to it as Mixed ChPT (MChPT). It is a 
kind of Partially Quenched theory which takes into account also the taste breaking of the 
staggered sea quarks. The valence quarks x — 1,2 and the ghost valence quarks x = 1,2 
have degenerate mass m m z, while the staggered sea quarks 4ii,4<i, 4s have masses m u = 
and m s . There is no unique recipe how to match the valence quark mass m va i and the sea 
quark mass m M , so the mixed theory always retains some features of partial quenching away 
from the continuum limit. 

The bubble contribution is obtained by applying the Wick contractions to © 



TjMChPT _ td2 

-°4+4+4 ~~ -°o 



^1 sea \^sea l) {^sea 2 1 2 sea) 
35=1,2,1,2 sea=4u,4dAs 

(13) 

To reduce four tastes per sea quark to one taste, the last term will be multiplied by 1/4 and 
the appropriate expressions for disconnected propagators from ^H] wm be used. 
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Figure 5: The effective mass of the scalar correlator for MILC simulation with 2 + 1 staggered sea 
quarks and staggered valence quarks. The effective mass is multiplied by r\ — 1.6 GeV -1 as in [3] 
for easier comparison. The dependence on am U/ /^ is shown, while am s is fixed to am s = 0.05 and 
am s = 0.031 on coarse and fine lattices, respectively. Triangles present effective mass of the bubble 
contribution in 2 + 1 Staggered ChPT Q at t = 20 for coarse MILC lattice and at t = 30 for 
fine MILC lattice; the corresponding error-bar denotes the variation due to the uncertainty in the 
hairpin parameters (|TU|) . In the continuum limit i? 2 5 +i l " PT © reduces to B ChPT and its effective 
mass is given by circles. Squares represent effective mass from MILC simulation jSJ. The relevant 



masses 2M 7T5 and M n5 + M V5 from MILC simulations are also shown (M, 



I/If 2 + 



The propagators for valence-valence mesons with x, x' = 1, 2, 1, 2 and mo — > oo are JB] 
= , 2 - 6 " 2 , M* aliVal = 2B m val , ei, 2 = 1 , e i 2 = -1 (14) 

K IV1 val,val 

1 (fc 2 + A^ J )(fc 2 + Aig J ) M^^M^+^A, 

<*» l *- v > dfa ' 3 (F + ML„J 2 (^ 2 + |M 2 7 + fMlJ ' Ml = Ml + a 2 A 7 ' 

while the propagators for valence-sea mesons with sea — u,d,s and x — 1, 2 are 

($ z seal^seaz) = , 2 ^ , ^lal, S ea = B (m va i + 17l sea ) + a 2 A Mix . (15) 

fc + M val,sea 

The propagators ()14I15|) depend only on two taste breaking parameters a 2 Aj and a 2 AMix- 
The taste-singlet breaking a 2 A/ of the sea-sea pion mass is taken from [5], while the taste 
breaking a 2 AMix of the valence-sea pion mass has not been determined yet and it is the only 
free parameter in the present paper. 
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The bubble contribution is obtained by inserting the propagators (jl4H5J) to (fT3j) 

Bjg h PT {t) = F.T.[ B^ hPT (p) ] ff=5 with (16) 

(P + M^ + M^) 
2+ i W „2^ 3 {k + p)2 + Km {k2 + M 2 alval ) 2 k 2 + + |M| 7 

1 

+ 




(* + P) 2 + ML U k 2 + M^, (* + v? + k 2 + Ml 



val,s 



Our result agrees in appropriate limit 9 with the scalar correlator from Ref. [17] . which 
considers mixed quark actions but not staggered sea. In the continuum limit (a 2 Aj — > 0, 
A-Mix 0) the expression (fTTpj) reduces to B ChPT ffTTjl once the valence and the sea quark 



a 



masses are tuned as m va i = m u . Away from the continuum limit, the contributions e~ 2Mnt 
from the diagrams in Figs. andEt do not cancel and the correlator has the effective mass 
2M n at large t. 

Fig. |H1 represents the bubble contribution (|16j) when valence and sea quark masses are 
tuned by matching M va i tVa i = Mu 5 which is used by LHPC Collaboration [7]. The expression 
ffTHl) is evaluated for MILC coarse lattices using M Us and V = 20 3 x 32 from LHPC [3 
and a 2 A j from [Hj. The a 1 'Amix is the only unknown parameter and it is varied in the 
reasonable 10 range —0.02 < a 4 AMix < 0.08. The bubble contribution is sizable for small 
m u / d and we find it is negative for a^Amix ^ —0.01. This negativity can be attributed 
to the fact that all the sea-sea pions (except 7r 5 ) are heavier than the valence- valence pion 
when M va [ yVa i = Mu 5 and the correlator goes negative for the same reason as in Partially 
Quenched QCD [Hj. At small m n /^, the scalar correlator © is dominated by the bubble 
contribution and is therefore negative for a large range of values a 2 AMi X - At large m^, 
the scalar correlator is dominated by e~ ma0t for moderate t, so it is positive irrespective of 
the sign of the bubble contribution. A lattice study of point-point scalar correlator at small 
quark masses offers a way to determine the unknown parameter a^Auix via Eq. (fTTpj) . 

If the valence and sea quark masses are tuned by matching the masses of the valence- 
valence pion and the heaviest sea-sea pion (M vaiiVa i = M Ur ), the bubble contribution is 
relatively small and positive for a^Auix ^ 0.04 (see Fig. |7j). Although the bubble contribu- 
tion could be slightly negative for a^Auix ^ 0.04, we expect the lattice scalar correlator to 
be positive in this case, since the small negativity in the bubble contribution is most likely 
out-weighted by the positive contribution e~ ma0t in (jSJ). 11 

Various choices of tuning the valence and sea quark masses have therefore large effects on 
the scalar correlator, although all the choices are equivalent in the continuum limit. UKQCD 
Collaboration in fact proposed to tune the valence and sea quark masses by considering the 
sign of the scalar correlator [7], but the proposal has not been used in practice yet. 

9 The results agree in the limit t — > oo, a 2 /S.Mix — * and m u = mid = m s . 

10 The natural limit for the magnitude of taste breaking is |a 4 Aj\/M| < a 4 A/ ~ 0.08. The additional 
restriction in case of negative a 2 &Mix is the positivity of the mass M va i iSea , which requires a Ajfu > —0.03 
for the simulation with am u / ( j/ ar ns = 0.007/0.05. 

11 The rough size of the ao-exchange in the point-point scalar correlator can be estimated from C a o(t) = 
(l6B 2 f 2 /m aQ )e- m " ot where f a0 ~ 0.08(3) GeV and m aQ = 1.58(34) GeV. 
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Figure 6: The bubble contribution B^{ hPT (t) (|TH|) for the simulation with chiral fermions on 
2 + 1 MILC staggered configurations with volume 20 3 x 32 and coarse a. The valence and sea quark 
masses are tuned by matching M va i jVa i = M\j h . The unknown parameter a^/^Mix is varied in the 
reasonable range -0.02 < a 4 A Mi x < 0.08. 



3 Conclusions 

The Nature of the observed scalar resonance ao(980) with / = 1 is not revealed yet. The 
determination of ao mass on the lattice could indicate whether ao(980) is a conventional qq 
state or perhaps something more exotic. 

The ao mass is conventionally determined from the flavor non-singlet scalar correlator 
on the lattice. However, the interesting contribution e~ ma0t to the scalar correlator is ac- 
companied by the contribution from multi-hadron intermediate states. The most important 
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Figure 7: The bubble contribution B^^ hPT (t) Q16|) for simulation with chiral fermions on 2 + 1 
MILC staggered coarse lattice with am u /,i/am s = 0.01/0.05 and volume 20 3 x 32. The valence and 
sea quark masses are tuned by matching M va i jVa i = M\j I . 



multi-hadron intermediate state is the bubble contribution B, which is the intermediate state 
with two pseudoscalars. In this paper we determined the size of the bubble contribution for 
simulations with three sea quarks of mass m u = rrid ^ m s . This offers a way to determine m a Q 
by fitting the lattice scalar correlator to C(t) = Ae~ ma0t + B(t). 

We determined the bubble contribution for the case of simulations with staggered sea and 
staggered valence quarks, as well as for simulations with staggered sea and chiral valence 
quarks. We found that the bubble contribution dominates the scalar correlator for small 
masses of u/d quarks. In proper QCD, the lightest two-pseudoscalar intermediate state is 
iTT] and the bubble contribution drops as e -( M ^+ M n) t a t large t. The effects of taste breaking 
or mixed quark actions make possible also the intermediate state with mass 2M n , so the 
correlator drops as e~ 2MlTt at large t. This unphysical contribution to the scalar correlator 
vanishes only in the continuum limit. 

• We used Staggered ChPT for predicting the bubble contribution in case of staggered sea 
and staggered valence quarks. Its sign depends on the values of the input parameters, 
indicating the broken unitarity of the theory due to the \/ Det trick. All the input 
parameters have already been determined for the case of MILC configurations and 
they render positive bubble contribution. 

• The case of staggered sea and chiral valence quarks was explored using Mixed ChPT. 
The size and the sign of the bubble contribution depends on one unknown parameter 
a 2 'Amzx- The bubble contribution and the scalar correlator are negative (positive) for 
most of the values —0.02 < a 4 A Mix < 0.08 if the valence and the sea quark masses 



13 



are tuned by matching the masses of valence-valence pion and the lightest (heaviest) 
sea-sea pion. The comparison of the point-point scalar correlator from lattice to our 
prediction offers a way to determine the unknown parameter a 2 A Mix- 

Our analytic prediction for the bubble contribution B(t) will be helpful for the determination 
of rriao by fitting the lattice scalar correlator to C(t) = Ae~ ma0t + B(t). Such extraction of 
scalar meson mass is possible for moderate quark masses and times where e~ ma(,t is not 
negligible with respect to the bubble contribution. 

Note added: 

Some of the results given in the present paper were already presented in the conference 
proceedings by the same author. 
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Appendix A: Coupling of a scalar current to two pseudoscalars 

Here we derive the coupling of a point scalar current d{x)u{x) to a pair of pseudoscalar 
fields at the lowest order of ChPT. The effective scalar current can be determined from the 
dependence of the QCD Lagrangian and the Chiral Lagrangian on the spurion field Aidu, 
where M. is quark mass matrix. The quark current d(x)u{x) is given by —dC.QCDldM.du-, so 
the effective current is represented by —dCchPr/ dAidu- The Chiral Lagrangian is 

CchPT = If^Ud^] + \B Q fTr[M^U + rfM\ + ... , (17) 

where U = exp{\/2i§ / f] incorporates the SU(3)l x SU(3)r pseudoscalar field matrix $, 
/ ~ 95 MeV, while the dots indicate terms of higher order or terms which are independent 
of M.du- The effective current is therefore equal to [TT] 

d{x)u{x) ~ -^B f[U(x) + U\x)] ud = S [$ (x) 2 } ud + - , (18) 

so the coupling of the current d(x)u(x) to a pair of pseudoscalar fields $ is equal to the slope 
parameter Bo = M 2 /(2m q ). The result (|THj) is derived within conventional ChPT here, but 
it applies also to various extensions of CHPT (Staggered, Mixed, Quenched and Partially 
Quenched) at the lowest order. In the extended case $ in Eq. (JTSj) represents the field matrix 
of the extended effective field theory. 



Appendix B: Pseudoscalar propagators in Staggered ChPT 

Here we list the propagators of pseudoscalar mesons in SChPT [Tlj for three sea-quark 
flavors (with one or four tastes) in case of m u = rrid ^ m s . All the propagators in Euclidean 
space-time have the connected part 

(4> b ff'\<t> b f'f)con = fc2+ 1 M 2 ; > ( 19 ) 

where Mff is the mass of the un-mixed meson with flavor //' and taste b = 1, .., 16 

M 2 ff , b = B (m f + m r ) + a 2 A b , M Ub = M Wb = M UUb , M Sb = M SSb , A 5 = . (20) 

The taste- vector, taste- axial and taste-singlet mesons have also the disconnected part. The 
disconnected parts of (4> b uu \<f> b u1l ) , ($a|$tf) and {<f> b uu \4>dd) d are all equal to 



i i 4 k 2 + M 2 Si 

(<Puu\<t>uu)d lS c = ~- (p - M 2 j)(fc 2 + 1^ + 2 M | } fOTm 



OO 



h 2 4- M 2 

lV \±V v _ 2 5 , fc +M Sv 

\uWuu)d,sc ~ V{k2 + M ^ )(p + M ^ )(p + 

^2 _|_ M 2 

^tUsc = -^^A^ ~ M 2J {k 2 + M2J {k 2 + M 2J ■ (21) 
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These hairpin couplings mix the flavor neutral mesons with /, V, A tastes and the resulting 
mass eigenstates are [13] 



M ni = M Vl , Ml = ±M£ 7 + \M 2 Sl , M v ,=oo for m oo 

M^ VA = M UvA , M WA = m v (C a 2 5' VA ) , M <A = rry(C a 2 5' VA ) (22) 

with C = 1 in four-taste theory 2-4 + 4 and C = 1/4 in a single-taste theory 2 + 1 ^3] . The 
functions m v and m,/ are defined as 

171,(5) = ^l[M 2 v + Ml v +35-Z(5)} (23) 
ny(cT) = yJ\[M 2 v + M| v +35 + Z(<5)] with 
Z(<5) = ^/(M| v - M£~P - 25(M| v - M 2 v ) + 9^ 

for the vector taste and analogously for the axial taste. 



Appendix C: The sign of bubble contribution in SChPT 

Here we explore the sign of the bubble contribution B SChPT (t) given in (jSJ). This is relevant 
since the positive definite bubble contribution to scalar correlator is expected in a proper 
unitary quantum field theory We consider the unitary four-taste theory 2-4 + 4 and a 
single-taste theory 2 + 1, where the later loses unitarity due to the fourth- root trick. 

Since the sign depends mostly on the values of the hairpin couplings a 2 5' v and a 2 5' A , we 
concentrate on a special case with m u = = m s , Ai 
simplicity. In this case 



and a 2 5[, = a 2 5' A = a 2 5' for 



lim B SChPT (t) 



Bl 



f — >oo 



2ttL 3 



4 N 3tt 



-2M^t 



8 a 2 5' f dis (C a 2 5" 



fdis{5) 



7T 



6 5 M 2 Mjf 



M v ,e 



-2Mrt 



(24) 



M v , = ^/Ml + 35 



with C — 1 in four-taste theory and C = 1/4 in a single-taste theory. It is easy to find 
a positive value of a 2 5', which gives B2^ PT (t) < indicating the breakdown of unitarity 
due to the fourth root trick. On the other hand, we find that 5|^f T (t) > for any values 
of a 2 5', Mn and t. The positive definiteness of B^+^if) can be demonstrated by noting 
that 5!£h T (0) > 0, BfS+f^oo) = and dB$$^ T /dt < for all t, which prohibits any 
extremum with a negative value. The positive definiteness of the scalar correlator in the 
four-taste theory is a consequence of the unitarity. 
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